Reduced Joule heating in nanowires 
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The temperature distribution in nanowires due to Joule heating is studied analytically using a 
continuum model and a Green's function approach. We show that the temperatures reached in 
nanowires can be much lower than that predicted by bulk models of Joule heating, due to heat loss 
at the nanowire surface that is important at nanoscopic dimensions, even when the thermal conduc- 
tivity of the environment is relatively low. In additition, we find that the maximum temperature in 
the nanowire scales weakly with length, in contrast to the bulk system. A simple criterion is pre- 
sented to assess the importance of these effects. The results have implications for the experimental 
measurements of nanowire thermal properties, for thermoelectric applications, and for controlling 
thermal effects in nanowire electronic devices. 



The thermal properties of nanostructures such as 
nanowires, nanotubes, graphene, and nanometric devices 
have attracted much attention recently because of the 
unexplored physics of thermal transport at nanoscale 
dimensions. In particular, the thermal properties of 
nanowires (NWs) have been studied both for possi- 
ble applications and for basic science, and an impor- 
tant question is that of the thermal properties during 
current flow. Indeed, Joule heating leads to reduced 
performance P, and breakdown [!, [j] of nanowire elec- 
tronic devices. It is also utilized to probe NW electronic 
transport properties @, and to controllably functionalize 
nanowires [fj|. Joule heating is also relevant to applica- 
tions in thermoelectrics[7H9(, to current- heated magnetic 
nanowires [10| . and is important in assessing the reliabil- 
ity of field-emission devices based on NWs. In all of these 
cases, the temperature distribution along the NW due to 
Joule heating plays an important role, and so far analyt- 
ical expressions for the full nanowire geometry in contact 
with a thermal environment have not been presented (Tl). 

Here we provide such expressions and show that the 
nanometer size plays a crucial role in determining the 
temperature distribution in Joule-heated nanowires, due 
to heat flow from the surface that cannot be neglected 
due to the high surface to volume ratio, even when the 
thermal conductivity of the environment is relatively low. 
As a consequence, the temperatures reached inside the 
nanowire are significantly lower than predicted by bulk 
models, which are often used to model temperature dis- 
tributions in NWs. The maximum temperature also 
scales differently with the length of the NW, increasing 
much more weakly with length. A criterion is proposed 
to assess the importance of these effects, which can be 
used, for example, when evaluating the need for more in- 
volved three-dimensional temperature simulation tools, 
or when designing thermal measurement systems. 

To be specific, we consider the system shown in Fig. 
la. There, a nanowire of radius R and length L is ther- 
mally and electrically contacted by two infinite planar 
contacts at z = ±L/2 and held at temperature T = 0. A 



uniform electronic current I flows through the NW, caus- 
ing Joule heating. In the steady-state, the temperature 
distribution between the two contacts is given by 



d 2 T (r, z) 1 dT (r, z) d 2 T (r, z) 
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where p — I 2 p/(itR 2 ) 2 knw with p the NW electrical re- 
sistivity and knw the NW thermal conductivity. The 
key physics is found in the presence of the theta func- 
tion on the right hand side of this equation, and in the 
boundary conditions at the surface of the NW. These 
boundary conditions can take different forms depending 
on the thermal relation between the nanowire and its 
environment [l2l| . Here we focus on the case of intimate 
contacts, but we expect the results to be relevant to a 
broad range of situations. In addition, the case of in- 
timate contacts has the advantage of not requiring un- 
known interface parameters, and is relevant to core/shell 
NWs with minor modifications to account for the shell 
thickness. 

For intimate thermal contacts, the temperature is con- 
tinuous at the interface Tnw {?, z )\ t= r — T env (r, z)\ r=R 
, and because the environment surrounding the NW has 
in general a different thermal conductivity K env the flux 
continuity at the surface gives a boundary condition 



knw 



dT NW (r, z) 



dr 



dT env (r, z) 



r=R 



dr 



(2) 



r=R 



When K env = 0, the temperature gradients are entirely 
in the axial direction and the problem reduces to the 
steady-state temperature distribution in a bulk material 



d 2 T bu i k (z) 



-p. 



(3) 



This equation is often referred to as that for "one- 
dimensional" heat transport, and is often used to analyze 
the thermal properties of NWs. The solution of the bulk 
equation for a NW between z = —L/2 and z = L/2 is 
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FIG. 1: Panel (a) shows the system under consideration: a 
nanowire of length L and radius 7? between two plates held 
at T = 0. The nanowire and environment have thermal con- 
ductivities knw and n env . The color map shows a typical 
temperature profile with blue corresponding to zero tempera- 
ture and red to the largest temperature. Panel (b) shows the 
temperature distribution along the length of a 10 nm radius 
nanowire, plotted for three values of the ratio of thermal con- 
ductivities between the environment and the nanowire. Panel 
(c) is the radial temperature distribution normalized to the 
temperature at the origin, for the same nanowire as in (b). 



To obtain the solution for the NW, we first calculated 
the Green's function for the differential equation in Eq. 
([T]) with the boundary condition in and integrated 
over the source term p6 (r — R). The solution is readily 
obtained as 



T(r, z) 
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where f3 — n env / 'knw ', and J„ and K v are modified 
Bessel functions of order v. In the limit (3 = the solu- 
tion can be shown to reduce to that of the bulk case, Eq. 
©• 

Figure la shows the calculated temperature in the 
whole system for a general case, showing the heat loss to 
the environment, while Fig. lb shows the temperature 
profile along the axis of a 10 nm radius NW of length 1 
/xm for three values of the ratio of thermal conductivi- 
ties j3. While for small values of /3 the profile approaches 
that of the bulk solution, for larger values of f3 the maxi- 
mum temperature reached in the nanowire can be orders 
of magnitude lower. (This arises despite the fact that 
the radial temperature gradients in the nanowire are still 
fairly small, as shown in Fig. lc.) Figure 2 shows the 
temperature reduction more clearly by plotting the max- 
imum temperature as a function of /3 for NWs of length 
1 fim and of radii 1, 10, and 100 nm. The deviation from 
the bulk solution depends strongly on the NW diame- 
ter, with the maximum temperature dropping rapidly as 
Kenv/nNW increases beyond a value that is strongly di- 
ameter dependent. 

As an example, for a 10 nm radius NW with the ther- 
mal conductivity of bulk SiGe (10 W/mK) in air with 
thermal conductivity 0.045 W/mK, we have /? = 0.0045 
and the maximum temperature is almost a factor of five 
smaller than the bulk solution would predict. The same 
NW surrounded by Si02 would have /3 ~ 0.1 and the 
effect is even more striking. (While the actual thermal 
boundary conditions for these systems might be differ- 
ent than that given by Eq. (2), the same behavior is 
anticipated with (3 replaced by the equivalent parameter 
describing the thermal properties of the interfaces). 

As expected, these effects become even more impor- 
tant for longer nanowires, with a further reduction of the 
nanowire maximum temperature. However, in addition 
to this qualitative expectation, there is also a change in 
the scaling of the maximum temperature with length. In- 
deed, as Fig. 3 shows, the bulk solution predicts a maxi- 
mum temperature that scales with the length squared; in 
contrast, for a nanowire with larger /3, the scaling is much 
different, and depends much more weakly on length. 

To provide a simple expression detailing these effects 
we consider the first mode in the above expansion, nor- 
malized by the first mode of the bulk expression: 
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FIG. 2: Maximum temperature in Joule-heated nanowires as 
a function of the ratio of the thermal conductivities of the 
environment and the nanowire. The dashed lines are the pre- 
dictions from Eq. (8). 



Expanding for small irR/L we obtain for the temperature 
in the middle of the nanowire 



T(0,z) 
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valid for /3 > 



— i | In (2^5) | . This expression is plotted 
in Fig. 2 for the NWs of radii 1 and 10 nm; while the mag- 
nitude of this expression is higher than the actual solu- 
tion, it nevertheless represents the functional dependence 
on ft very well, i.e. it decreases as 1//3. Furthermore, the 
dependence on L is also well represented by this expres- 
sion, as the dashed line for the R — 10 nm NW of Fig. 3 
shows. By including the L 2 dependence of the bulk so- 
lution, we get from Eq. © the scaling T"j^ ~ |ln (^) | 
which shows the much weaker dependence on L. From 
this analysis we deduce that the thermal conductivity of 
the environment at which the bulk solution ceases to be 
applicable is 



Kenv £ K NW 



(9) 



indicating the importance of the aspect ratio R/L and 



thus the nanoscale dimensions. For a 10 nm radius NW 
of length 1 /im, deviations arise even when n env is only 
0.1% of K NW . 

In summary, we presented a simple analytical expres- 
sion for the temperature distribution in Joule-heated 
nanowires. We find that the temperature reached in 
nanowires during Joule heating can be significantly less 
than that obtained from standard bulk models. This 
originates from heat loss at the nanowire surface that 
cannot be ignored due to the high surface to volume ratio, 
even when the thermal conductivity of the environment 
is relatively low. The work presented here is straightfor- 
ward, but given the ubiquitous presence of Joule heating 
in nanostruc tures, it should be relevant to many a reas 
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FIG. 3: Maximum temperature in Joule-heated nanowires as 
a function of the length of the nanowire. The temperature 
is normalized by the maximum of the bulk solution, which 
scales as L 2 . 



of nanoscale thermal transport, including thermoelectric 
devices based on nanowire arrays, nanoscale transistors, 
and electrically driven light-emitting nanodevices. 
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